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ABSTRACT 

We present results on the spin-independent quark- antiquark potential in SU(3) 
gauge theory from a simulation on a 48 3 x 64 lattice at j3 = 6.8, corresponding to 
a volume of (1.7 fm) 3 . Moreover, a comprehensive analysis of spin- and velocity- 
dependent potentials is carried out for SU(2) gauge theory, with emphasis on 
the short range structure, on lattices with resolutions ranging from .02 fm to 
.04 fm. 
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1. Introduction 



Quarkonia spectroscopy provides a wealth of information and thus constitutes an 
important observational window to the phenomenology of strong interactions. Within 
their non- (or semi-) relativistic setting, potential models have been proven to describe 
the empirical charmonium and bottomonium spectra remarkably well. The resulting 
phenomenological potentials can be compared to lattice predictions from the QCD 
Lagrangian. 

Lattice QCD techniques offer on the other hand a direct access to the hadron 
spectrum, without recourse to the potential picture. In the heavy quark sector (with 
quark mass M), however, the full fledged lattice calculation would require (at present) 
prohibitively fine lattice resolutions a < M -1 . It is therefore more practical to expand 
the QCD Lagrangian in powers of M~ l into an effective non-relativistic QCD action 
(NRQCD) which can be evaluated subsequently by use of lattice methods.^ 

The traditional semi-relativistic potential model approach to spectrum calcula- 
tions assumes the validity of the instantaneous approximation. This approximation 
can be established from first principles lattice computations by comparing spectra 
and wave functions as obtained from NRQCD to results from lattice potentials. 

Pioneering attempts to compute corrections to the static potential on the lattice 
have been launched in the mid eighties.@i In the meantime, tremendous progress has 
been achieved, both in computational power and methods. In view of the general 
interest in the potential formulation of quark binding problems, it appears to be 
timely for a new lattice determination of spin-dependent (sd) and velocity-dependent 
(vd) potentials. 

As a first step into this direction, we present recent results on the static quark 
antiquark potential from SU(3) gauge theory, as well as results on sd and vd potentials 
from SU(2) gauge theory. The limitation to two colours will not yet allow to proceed 
to spectrum calculations but will hopefully disclose the key features of confinement 
at work. 

2. The Central Quark- Antiquark Potential 

In the limit of infinitely heavy colour sources (M — > oo), the Born-Oppenheimer 
approximation can be applied and — after integrating out the gauge degrees of free- 
dom — the underlying QCD binding problem becomes nonrelativistic. 

Theoretically, one expects the leading order (in M~ l ) potential to be dominated 
by one-gluon exchange at short distances, Vo(r) = —Cf^^-, with a running coupling 
a(r) which depends logarithmically on r. On the other hand, for large separations 
r, one would expect a bosonic string model to hold which predicts the asymptotic 
behaviour Vo(r) = kt — + • • • (with string tension k). A priori it is not clear 
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Fig. 1. The static SU(3) QQ potential, obtained on a 48 3 x 64 lattice at = 6.8, corresponding to 
a lattice spacing a « .035 fm. 



when large or small r behaviour would set in and which form the potential takes in 
the region of intermediate r. 

Phenomenologically, the Cornell potential V^ eri (r) = nr — - has been proven to 
describe the spin averaged charmonium and bottomonium spectra successfully^ with 
one set of parameters, k and e, in accord with the flavour independence of strong 
interactions. It is well known, though, that the binding energies of low lying states 
are not overly sensitive to the shape of the potential at distances smaller than 0.2 fm 
or larger than 1 fm. 

Therefore, a lattice determination of the QCD potential in the intermediate r 
region is of particular interest. In Fig. [L] we show high statistics results on the central 
potential from SU(3) gauge theory on a 48 3 x 64 lattice at (3 = 6.8, which corresponds 
to a lattice resolution a « 0.035 fm. The data for r > 0.3 fm has been fitted to a 
Cornell type parametrization (solid curve). At small r we observe deviations from 
this parametrization. This is in accord with running coupling effects.il Note, that in 
the lattice computation vacuum polarization due to sea quarks has been neglected 
and thus no string breaking occurs. 

Before one can proceed to predict spin-averaged quarkonia spectra from such 
lattice potentials, one has to investigate possible corrections to the infinite mass limit 
at realistic charm and bottom quark masses. In solving the Schrodinger equation with 
our potential, we find the average speed of the quarks (in units of the speed of light) 
to be (v 2 ) ~ 0.25 and (v 2 ) ~ 0.09 for the charmonium and bottomonium ground 
states, respectively. This leads us to expect that, at least in the case of charmonium, 



the phenomenological potential, lp/ ie „(r), (which has been tuned to reproduce the 
empirical spectrum) might deviate by substantial 0(v 2 ) corrections from the static 
potential, Vo(r). Needless to say, that fine and hyperfine splittings have their origin 
in such v or, equivalently, 1/M corrections. 



3. 1/Af-Corrections 

Starting from a Foldy-Wouthuysen transformation of the Euclidean quark prop- 
agator in an external gauge field, the asymptotic {t — > oo) expression W(r,t) oc 
exp(— Vo(r)t) can be derived for the static potential, V (r), where W(r,t) denotes 
the expectation value of the familiar Wilson loop with spatial extent r and temporal 
extent t. Perturbing the propagator in terms of the inverse quark masses Mf 1 and 
M^ 1 around its static solution, one arrives at the semi-relativistic Hamiltonian (in 
the CM system, p x = — p 2 = p), 

P 2 / l n (p 2 ) 2 ( l l 



i/ =tU + ^J-^(w + wJ +l/(r ' L ' Sl ' S2 ' p) • (1) 

where the potential consists of a central part, sd and vd corrections, 

V(r,L,S 1 ,S 2 ,p)=V (r)+V sd (r,L,S 1 ,S 2 )+V vd (r,p) . (2) 

The sd contributions have been derived by Eichten, Feinberg and Gromesi while the 
vd terms have been elaborated by Barchielli et a/.:Q 

vWr,L, Sl ,s,) = h^^ W{r) + 2V ; {r)) + h^i vi{r) 

+ f^^-^W + ^M (3) 
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with Si = Si /Mi, 1 ; = Li I Mi and 



VUr.p) = l(^ + ^ 2 (Mr) + V a (r)) (4) 
+ L\, v iL jVi (r) + ( S 4- r ^\v c (r)\) (5) 
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+ E {4, v{, (5 i3 V d {r) + - 7?) }^ (6) 

with Vi = Pi/Mi, respectively. {-,-,-}w denotes Weyl ordering of the three arguments. 

The sd and vd potentials V[, V^', V 3 , V4 and V a - V e can be computed from 
lattice correlation functions in Euclidean time, C(T). This is done by measuring 
expectation values of Wilson loops with two colour field insertions (ears) within the 



temporal transporters, divided by the corresponding loop without ears (T denotes 
the ear to ear temporal distance). One ear excites the gluon field between the two 
charges while the second one returns the field into its ground state. To obtain the 
potentials, an integration over all possible interaction times (temporal positions of 
the ears) has to be performed. 

The minimal distance of an ear to an "end" of the Wilson loop, occurring within 
the integration, AT, represents the time the gluon field has at its disposal to decay 
into the ground state after creation and, therefore, governs excited state contamina- 
tions. The spatial transporters within the Wilson loops are smearedi to suppress such 
pollutions from the beginning, allowing us to reduce AT. As an additional technical 
trick we integrate out temporal links analytically^ which results in reduced statisti- 
cal fluctuations. Finally, we exploit transfer matrix techniques to obtain asymptotic 
results from a finite integration range (in T). Details about this procedure will be 
published elsewhere. 

4. Matching of the Effective Hamiltonian to QCD 

The sd and vd potentials are computed from amplitudes of correlation functions 
rather than from eigenvalues of the transfer matrix. This gives rise to renormalizations 
in respect to the corresponding continuum potentials. A different way to illustrate 
the necessity of renormalization is the fact that the colour electric (magnetic) "ears" , 

9F,u = ^ 2 {U, u -Ul)+0{a 2 ) , (7) 

which have been inserted into Wilson loops, explicitly depend on the lattice scale, a, 
and discretization. 

However, renormalization is not a pure lattice problem in this case. By truncating 
the expansion of the QCD Lagrangian in powers of M -1 at a given order, the ultra- 
violet behaviour is changed in respect to the full theory. Therefore, the effective 
Lagrangian has to be matched to full QCD at a renormalization scale /x < M" giving 
rise to renormalization constants q(/x, M), connecting a QCD potential Vi(r;M) to 
the corresponding potential, computed in the framework of the effective theory, e.g. 
Vi(r; M) = Cj(/i, M)Vi(r; /x, M). This problem, which becomes visible beyond the tree- 
level, has been approached systematically for sd potentials in the context of heavy 
quark effective theory by Chen et aL0 

Up to an additional <5 3 (r)-like factor, that originates from mixing of the dimension 
6 spin-spin interaction term with a four fermion contact term, a result is obtained 
which can be rewritten into Eq. (^) by substituting the naive potentials V{ by renor- 
malized potentials Vi ren which can be related to each other by a renormalization 



a On the lattice, /x ~ w/a. 




Fig. 2. Test of the HM renormalization procedure against the Gromes relation. The data on V% — V{ 
correspond to results from two different lattice spacings. The fit curve refers to the central force V$ 
as computed from smeared Wilson loops. 



constant, c 3 b , 

V{, ren (r) = c 3 (/i,M)V 1 / (r, / i) + (c 3 ( / i,M)-l)V '(r) , (8) 

V^rJr) = <%(/*, M)V?(r.M) , (9) 
V 3/i , ren (r) = c 2 3 (fi,M)V 3/4 (r,fi) . (10) 

The constant C3 is known to one-loop perturbation theory.0 As C3 is likely to 
be dominated by higher order perturbative and nonperturbative uncertainties in the 
low energy regime of interest, we apply the (nonperturbative) HM renormalization 
procedure induced by Huntley and MichaelS This procedure has originally been 
invented to remove the continuum-lattice renormalization problem but — as a by- 
product — also cures the matching problem of the effective heavy quark theory. The 
success of this approach can be checked numerically in two ways, namely (a) by 
varying the lattice resolution a and checking scaling of the results and (b) by testing 
against the Gromes relation^ between spin-orbit potentials and the central potential 
(which does not undergo renormalization)^ 

K(r) = V'(r) - V'(r) . (11) 



fc For simplicity we assume M — Mi = M2. m 
c Two additional constraints, relating vd potentials to the central potential, have been foundLl: 
rV^{r) = -2V c {r) - AV e {r) and rV '(r) - 3Vo(r) = 6V b {r) + 12V d {r). 
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Fig. 3. Scaling plot of the long range spin-orbit potential V{ (in units of the string tension, n). The 
curve corresponds to a one-parameter fit of the form V{{r) = — n — a/r 2 . 
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Fig. 4. The spin-spin potential V4 in lattice units. 
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Fig. 5. The velocity-dependent potential Vd in lattice units at j3 = 2.74. 

The present simulations have been performed on 32 4 and 48 4 lattices at (3 = 2.74 
and (3 = 2.96 which corresponds to lattice spacings a ~ 0.041 fm and a ~ 0.020 fm, 
respectively. The above physical scales have been adjusted such that the string tension 
comes out to be a/k = 440 MeV. In Fig. |^ we check our data on V 2 ' — V{ (in units of 
the string tension, k) against the force, obtained from a fit to the central potential, 
Vo(r). As can be seen, the two data sets scale nicely onto each other and reproduce 
the central potential up to lattice artefacts at small r. 

5. Results 

In the continuum, tree-level perturbation theory yields the following expectations 
for the spin-orbit and spin-spin potentials, 

P tP 

Vl = , V^ = - , V 3 = - , V A {r) = 8n5 3 (r) , (12) 

where e = C^a and CV = 3/4 for SU(2). Note, that the first spin-orbit potential does 
not contain vector exchange contributions and, due to the Gromes relation, should 
be of the form (assuming Eq. (0) to be valid and a Cornell parametrization of the 
central potential), V[ = —k. 

The data for V{ is displayed in Fig. [| Our lattice resolution enables us to establish 
an attractive short range contribution to V[ that can be well fitted to a Coulomb (1/r 2 ) 
form (in addition to the constant long-range term, which is in agreement with the 
string tension, k). This term amounts to about one quarter of the Coulomb-like 




contribution to the static potential. 

We confirm the second spin-orbit potential V 2 ' to be definitely of short range 
nature. Moreover, up to the first data point, it qualitatively agrees with tree-level 
lattice perturbation theory. This agreement can be made quantitative by allowing for 
a running coupling parameter. The same holds true for the spin-spin potential V3. 

The remaining spin-spin potential, V4, exhibits oscillatory behaviour as a lat- 
tice artefact (Fig. |j) and can largely be understood as a ^-contribution, according 
to Eq. fljjP . The tree- level lattice perturbative expectation 87rc5£(r/a) is indicated 
by squares in the figure. The normalization has been obtained from a c/r 2 fit to 
the V 2 ' data points. The error bands without symbols are obtained by using lattice 
single gluon exchange with an infra-red protected two-loop running coupling in mo- 
mentum space.0 The range corresponds to different reasonable choices of the QCD 
A-parameter. Apart from the dominant £-like contribution, another very short ranged 
contribution seems to exist. 

Very recently, we have computed the vd potentials for the first time. We find 
reasonable signals for long range forces, as illustrated in Fig. [5]. 

6. Summary and Conclusions 

We have studied central and spin-dependent forces in SU(2) gauge theory in a 
high statistics lattice simulation. We find reliable renormalized potentials with good 
scaling behaviour. There is clear evidence for a short range scalar exchange contri- 
bution in the first spin-orbit potential at the level of 20-25 % of the Coulomb part of 
the central potential. The other sd potentials are found to be short ranged and are 
well approximated by perturbation theory. 

We are encouraged from the results of a feasibility study of velocity- dependent 
potentials. An extension of the present investigations to the case of interest, SU(3) 
gauge theory, is in progress. 
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